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Abstract 

Diffusion of ions through a fluctuating polymeric host is studied both by Monte Carlo 
f^. simulation of the complete system dynamics and by dynamic bond percolation (DBP) 

ON \ theory. Comparison of both methods suggests a multiscale-like approach for calculating 

' the diffusion coefficients of the ions. 
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1 Introduction 

Chain polymers carrying electro-negative atoms (e.g. oxygen or nitrogen) in their re- 
peat unit can act as solvents for certain salts. Well-known examples are Li-salts dissolved in 
polyethylene-oxide (PEO). At temperatures sufficiently above the glass-transition tempera- 
ture these polymer-salt solutions show significant DC ionic (Li + ) conductivities. Such "poly- 
mer electrolytes" offer widespread applications in batteries, sensors and fuel cells. From the 
scientific point of view, an important goal is to improve our understanding of the electrical 
conduction mechanism in polymer electrolytes, and in particular to elucidate the interplay 
between ion diffusion and the polymer network dynamics. 

Dynamic Monte Carlo (MC) simulation of the diffusion coefficient of few particles (ions) 
in a rearranging environment of polymer chains is hampered by the need to move every 
monomer (polymer bead or ion) with the same probability. Therefore most of the computa- 
tional time is spent moving the polymeric host without affecting the ionic configurations. On 
the other hand, earlier studies indicate that important features of the ion diffusion within a 
dynamical matrix of chain molecules can be described by a more coarse-grained model[|T]]. 
The idea is to map the diffusion process onto DBP-theory [0,0] and to determine the central 
quantity entering this theory, the renewal time distribution ip(t), from the time dependence 
of the local occupational correlation function due to the polymer chain dynamics. In this 
communication we perform tests of such a procedure for a hard-core lattice gas and for tracer 
diffusion in systems of athermal lattice polymers. It turns out that the DBP-concept compares 
favourably with simulations of the complete system dynamics and thereby saves about one 
order of magnitude in computer time. 



2 Hard-core lattice gas 



In order to explain our procedure we treat the simple case of a hard-core lattice gas, i.e. 



diffusion constant D(c) as a function of concentration c is known to a high degree of accuracy 
via dynamic pair approximations which yield a tracer correlation factor 



f( v D(c) l + (cos9) 

71 J ' Do(l-c) l-[(3-2c)/(2-c)]<cose)' 

with (cos 9) ~ —0.209 for a simple cubic lattice. 

Now we consider a lattice with static disorder in which randomly chosen sites are blocked 
and thus not accessible for the tracer particle. For this standard percolation problem we de- 
note the mean square displacement of the tracer by (r 2 (t)) . If on the other hand the blocked 
sites are globally and instantaneously rearranged according to an arbitrary waiting-time dis- 
tribution ip{t), generalized DBP-theory yields the following diffusivity at zero frequency in 
d = 3 dimensions [Q] 

_ l / °°d^(t)(r 2 (t))o 
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Our aim is to map the complete system of coupled ions and lattice chains onto this coarse- 
grained DBP-model. While (r 2 (t)) can be obtained in a straightforward manner from sim- 
ulations of ion diffusion in the frozen network, the determination of ip(t) requires more ex- 
planation. Within the spirit of previous work on the hard-core lattice gas we propose to 
determine ip(t) from the local occupational correlation function (nj(t)nj(O)) of a site i ad- 
jacent to a fixed tracer position. Let us introduce the probability $(£) that no renewal takes 
place within the time interval [0, t] after a previous renewal at an arbitrary t < 0. The joint 
probability (nj(t)ni(O)) that the lattice point i is occupied at t = and at t (not necessarily by 
the same particle) consists of two distinct contributions. The first one is the probability that i 
is occupied at t = and that no renewal occurs until t, which is given by c $(£). The second 
contribution describes the situation where one or more renewals have taken place until t. The 
corresponding probability is given by c 2 (l — $(£)). Hence $(£) is related to (rij(t)nj(O)) via 

c(l - c) 

According to 0], ip(t) = $"(t)A _1 where A = J 00 tip(t)dt denotes the mean renewal time. 
Equation @ thus can be rewritten as 

l f °°dt<t>"(t){r>(t)) 
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where we obtain $(t) from (|3|) by means of MC simulations for (nj(t)ni(O)). The resulting 
$(£) can be fitted with sufficient accuracy by a combination of P exponential functions, 

p 

$(t) = ^a ie xp(-A,t) (5) 
i=i 

On the other hand, (r 2 (t)) can by expressed in terms of the simulated mean square displace- 
ment (r 2 ) of a tracer particle making n steps in a frozen network via 

(r 2 (t)) = lira V expf-u^Kr^n (6) 
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Figure 1: Test of DBP-theory against simulations of the complete system and the dynamic 
pair approximation equation ([I]). The quantity shown is the tracer correlation function f(c) as 
a function of concentration. Also plotted are results obtained recently form a many-particle 
effective medium approximation [§]. 



where lu is the attempt frequency of the tracer. Inserting (Q) and (§) into (Q) and carrying out 
the integration one gets 

- = - lim T (^) 2 a, Y „ , (r "\ - (7) 
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Independent calculations for c — > 1 showed that in order to keep the error caused by the 
finiteness of iV and p below the statistical error, values of P = 3 and iV = 40 were sufficient. 
Note that the main decay of $(£) is roughly represented by a single exponential with decay 
rate A = —(d®/dt) t=0 = (5/Q)lu . 

In Figure (|l]) we compare the results of this theory with MC simulations (see also [0]) 
and the dynamic pair approximation equation (|T]) and also with an existing effective medium 
approximation (EMA) for a percolating lattice renewing itself with only one rate constant 
A = 5/6cjo [§p- As one can see, using more than one rate constant in the representation of 
$(£) according to (§) improves the results significantly. This observation suggests to apply 
the generalized DBP-theory to the situation described in the next section where the non- 
exponential character of the renewal process is even more pronounced. 

3 Results for athermal chains and conclusion 

In the case of a hard-core lattice gas the procedure described bears no computational 
advantage over existing methods but was merely considered as a test case. The situation 
changes, however, when we apply our approach to ion diffusion in a polymer network. No 
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Figure 2: Comparison of the present theory with full MC simulation using ion concentrations 
c; on = 10~ 3 (one ion) and c\ on = 10~ 2 (10 ions). For comparison the result of equation ([]]) for 
a hard-core lattice gas is also included. 



ulations of D become much more demanding because of the internal degrees of freedom of 
the host molecules and the larger statistical errors connected with the small concentration of 
tracer particles. 

Let us consider, for example, a system where the chains are made of ten beads, assigned to 
lattice sites, and linearly connected via nearest-neighbour bonds. Apart from site exclusion, 
which mimics a hard-core repulsion, no explicit interactions between beads are assumed. 
The size of the simulation box is L = 10 and periodic boundary conditions in all directions 
are employed. For the chain dynamics we use the standard algorithm as described in Ref. 
[EHUJh while point-particles individually perform nearest-neighbour hops. 

To apply our theory we again determine $(£) in analogy to the hard-core lattice gas, see 
(|]), which is now highly non-exponential, indicating the importance of temporal correlations 
in the associated renewal processes. A superposition of three exponential functions in (|5]), 
however, still gives a good description of $(£). Mean-square displacements (r^)o of ions 
within a frozen network of chains are determined using N = 200 steps, and the diffusion 
constant is calculated as before via equation ([7J). In Fig.(|2|) we compare our results for the 
tracer correlation factor as a function of the total concentration of occupied sites, c, with full 
MC simulations. Due to the connectivity of chains blocking effects at a given c on average are 
reduced and the correlation factor f(c) is larger than in the case of point particles, at least as 
long as c < 0.8. More important from the computational point of view is the fact that DBP- 
theory very well agrees with the full simulations, but saves about one order of magnitude 
in computer time. Mapping that problem onto a coarse-grained DBP-model thus appears to 
constitute an accurate and efficient method for investigating diffusion through a dynamic host 
of chain molecules. To support this conclusion further work is being carried out, in particular 
under varying chain lengths. 
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